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ABSTRACT 

The critical behavior of the semi-infinite Blume-Capel and Blume-Emery-Griffiths( 

models is investigated in the pair approximation of the Cluster Variation Method. 
Equations for bulk and surface order parameters and n.n. correlation functions 
are given, from which analytical expressions for the second order bulk and surface 
critical temperatures are derived. The phase diagrams of the Blume-Capel model 
are classified, and the existence of a surface first order transition is discussed. This 
transition is shown to be, under certain conditions, slightly reentrant, and the be- 
havior of the surface order parameters and correlation functions is given for such 
a case. The extension of our results to the Blume-Emery-Griffiths model is briefly 
discussed. 
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1. Introduction 

The spin-1 Ising models have been originally proposed to describe the 
critical behavior of magnetic systems[l], and have subsequently been applied 
to He 3 -He 4 mixtures[2] and to multicomponent fluids [3]. The most common 
and most extensively studied among these models are the Blume-Capel (BC) 
modelfl] and the Blume-Emery-Grifhths (BEG) model[2], the former being 
but a special case, with no biquadratic exchange interaction, of the latter. 

These models are not exactly solvable in more than one dimension (one ex- 
ception being the BEG model in a particular subspace of the phase space[4]), 
but they have been studied over infinite (i-dimensional lattices with a lot of 
different approximation techniques, and their phase diagrams are well-known. 

In the last few years some attention has been devoted to the study of 
these models over semi-infinite lattices, with different couplings at the sur- 
face and in the bulk. Benyoussef, Boccara and Saber [5] have investigated 
the BC model in the mean field approximation (MFA), giving a complete 
classification of the possible types of phase diagrams, while Benyoussef, Boc- 
cara and El Bouziani[6] have carried out the same investigation on the BEG 
model in a real-space renormalization group (RG) framework. A few other 
papers have been devoted to the application of the semi-infinite BEG model 
to the study of surface superfluidity in He 3 -He 4 mixtures, both with RG[7,8] 
and MFA [9], and a particular investigation of the BC model, in the effective 
field approximation (EFA), in the case of equal bulk and surface exchange 
interactions is due to Tamura[10]. 

As in the case of the spin- 1/2 Ising model, it is possible to have an ordered 
surface even when the bulk is disordered. When such a situation occurs, it 
is customary to speak of extraordinary (i.e. bulk) and surface transition, 
and the temperature of the extraordinary transition is lower than the surface 
one. Otherwise, bulk and surface disorder at the same temperature and 
with the same critical exponents and the transition is said to be ordinary. 
Finally, in the limiting case between the two situations above, one has a 
special transition, with bulk and surface disordering at the same temperature 
but with different critical exponents. 
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While these results are well-established, the literature shows a contro- 
versy about the order of these transitions. Some papers, where MFA[5,9] or 
RG[7,8] are employed, indicate clearly that the surface transition (and thus 
also the extraordinary one) can be either second or first order; on the con- 
trary, there are authors who exclude the possibility of a first order surface or 
extraordinary transition, both in a RG[6] and in an EFA[10] approximation 
scheme. 

In this paper we investigate the semi-infinite spin-1 Ising models in the 
pair approximation of the cluster variation method (CVM), which has been 
introduced by Kikuchi[ll] and has subsequently been reformulated by An[12], 
who put it into a much simpler form resorting to the Mobius inversion. 

In order to allow for the spatial variation of the local quantities which this 
method introduces, we divide the lattice into an infinite set of layers which are 
parallel to the free surface: since this procedure leads to an infinite number 
of equations, we make another approximation, which consists in treating the 
first layer below the surface as if it was bulk; this way we obtain only two 
finite sets of equations: one for the bulk, which is not coupled to the surface 
one, and one for the surface, which contains the bulk variables as parameters. 

The use of the CVM has two main advantages: it allows us to obtain 
the correlation functions in a quite straightforward way and to distinguish 
easily between second and first order transitions, the former being located 
analytically. Furthermore, we overtake the problem of dealing with a great 
number of variables and equations by considering a minimal set of indipen- 
dent functions (the order parameters and the correlation functions) of the 
elements of the reduced density matrices, and by developing moreover a pro- 
cedure which allows us to solve the equations for the correlation functions, 
leaving us with only two coupled equations for the order parameters. These 
equations, when solved numerically by iteration, give always thermodinam- 
ically stable solutions, i.e. minima of the free energy. From the equations 
for the order parameters we deduce an equation for the second order critical 
temperature, while first order transitions are determined numerically. 

The analysis, in the whole phase space of the BC model, shows two new 
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types of phase diagrams and a reentrant first order surface transition. For 
the latter, we give qualitative conditions about the region of the phase space 
where it occurs and an example of behavior of the order parameters and of 
the correlation functions. 

The paper is organized as follows: in Sec. 2 we determine the CVM 
free energy, the equations for the order parameters and the second order 
critical temperatures; in Sec. 3 we describe our main results: the classification 
of the phase diagrams for the BC model and the analysis of the reentrant 
phenomenon in the first order surface transition. In Sec. 4 we discuss the 
extension of our results to the BEG model and some conclusions are drawn 
in Sec. 5. 

2. The model and the approximation. 

The semi-infinite BEG model has hamiltonian 

ph beg = -JsYl s * s j - jBj2 SkSi+As Yl s * + Ab E s l 

(ij) (kl) i k 

{ij) {kl) 

where k, I are site labels, Si is the ^-component of a spin-1 operator at site 
i, denotes a sum over all nearest neighbors (n.n.) with both sites lying 

{ij) 

on the surface, denotes a sum over the remaining n.n., and (3 = (/csT) -1 

{kl) 

(with ks Boltzmann constant and T absolute temperature). Js and Jb 
(both positive, since we study the ferromagnetic case) are reduced exchange 
interactions, while As and Ab are reduced crystal fields and Ks and Kb are 
reduced biquadratic exchange interactions, respectively at the surface and in 
the bulk. 

As we mentioned in the introduction, we divide our lattice into layers, 
parallel to the surface, labeled by an integer n, n = 1 being the surface layer. 
For the sake of simplicity in the following we will consider a simple cubic 
lattice with a (100) free surface. 
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The CVM is based on an approximate expression of the entropy of the 
model as a sum of contributions by all the elements of a set M. of maximal 
clusters and all their subclusters. In the pair approximation, one chooses as 
M. the set of all n.n. pairs. 

Thus our entropy will be the sum of contributions of three different kinds: 
one due to one-site clusters on layer n, denoted by ai; one from two-site 

(n) 

clusters with both sites in layer n, denoted by a\, and one from two-site 
clusters with a site in layer n and another in layer n + 1, denoted by <J^\ 
All these contributions will be weighted appropriately yielding, for the total 
entropy a 



where N$ is the number of sites in a layer and 7 ranges over all the clusters 
above. The weights a 1 can be determined by means of An's equations[12] 
and, as well as the multiplicities N^/Ng, are lattice dependent. In our case 
we obtain 



aW = -4, c4 n) = -5, n > 2, N^/N s = 1, n > 1 
a (n) = a ( ? ) = 1; N^/N S = 2, N^/N S = 1, n>l. (3) 



Finally, a 1 is the entropy associated to the cluster 7 and is given by 



where p 1 is the reduced density matrix for cluster 7, which has to be deter- 
mined by minimization of the free energy. 

Since the model hamiltonian is diagonal the reduced density matrices 
turn out to be diagonal as well. Furthermore, they must obey the following 
constraints: 




(2) 
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/c B Tr(p T In p 7 ), 



(4) 



Trp 



= 1, 




e 



e 






e 



e 
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where pp Se = (Se\p p \6e) and so on, with 5 and e taking values +, 0, — , and 
it is assumed for p^ that the first index refer to the site in layer n and the 
second one to the site in layer n + 1. 

These constraints allow us to re-express the density matrices as functions 
of a reduced (with respect to the set of the elements of the matrices their- 
selves) set of order parameters and correlation functions; this is the first step 
in lowering the number of equations one must deal with when minimizing the 
free energy. Upon defining the order parameters (( ) denotes now thermal 
average and stands for a spin in the nth layer) 

y[ n) = (S^) and y M = (sl n) \ (6) 

and the n.n. two-points correlation functions 

y (n) = {S (n) S (n)^ y >{n) = ^(n) ^(n+1) ) . 

y (n) = {S (n) S W\ y /(») = ( <,(n) ^(n+l)^ y J<») = (S^ 

y (n) = {S (n) 2 S (n)* h y /(n) = ( <j(n) 2 <,(n + l) ^ 

(7) 

one obtains the following expressions for the elements of the density matrices: 

(n) _|_ (n) 

In) _ V2 ±Vl (n) _1_„(»). 

Ps ± — 2 ' ' s ~~ y2 ' 



p (n) 



+ 8eyt ] + (5 + e)y[ n) 



Jn) _ -, , (n) _ 2 (n) 



P 00 

„Xn) _ (n) , (n) (n) 
(n) _^ (n) 1/2 ^5 ± ^1 

Pp ±0 Pp ± O 



P Se 



(n) /(n) 1 (n) _,_ /(n) 

(») _ ^2 — 2/5 =»= 3/1 =F y 4 



P ±0 



(n+1) /(n) , (n+1) _,_ //(n) 
(n) j/2 -1/5 =F2/4 



P 0± 



» 



= i + vT ) -vt ) -y\ 



(n+l) 



(8) 



' V 00 

where now 5 and e take values +, — . 

We are now able to write down our approximate expression for the reduced 
free energy density / = j3{U — Tcr)/Ns; we obtain 

f = A s yi 1] - 2K S yW - 2J S 

oo 

+ J2 [Ab V? - y ( 3 n) - 2K B y { f 

n=2 



oo 

-J2[j B y} n) + K B y't y 

n=l 

oo 

+ E[ 2TV (p? ) lnpW)+'&(p? ) lnp5? 

n=l 

oo 

- 4Tr ( P W InpW) - 5 £ TV (p(»> lnp<»>) 



(9) 



n=2 



Requiring that / is a minimum with respect to all the y's yields an infinite 
set of coupled equations, which we report making use of the notation X^ = 
lnpi n ' ) , Ap™"* = lnp p n ' ) , X^ = lnp p ?\ and indices have the same meaning as 
in (8): 



4(A p n) +0 - A p n) _ )+ 

/,(n) x (n) \( n - 1 ) 

l A P ' +0 _ V -o + V 



0+ 



A ?" 1) o-) + ^ n) ( A ^ ) +- A ^ ) -) = ° 



2-D n + 4(Ap _|_ + Ap 

(n) 



o\( n ) \ , 
ZA P 00/ " t ~ 



(V +., 



+ a!? } 



ai n) (Ai n) + + Ar j --2Ar )=0 



» 



-2Ai? ) oo + Al?- ,i 

(n) 



+ a!?- 1} 



-4J n + (A p n) ++ + Ap n) __ - 2Ap n) + _) = 

-4J B + (A<?> + X^__ - X<£\_ - X^_) = 



Ap 



(n) 



- Ap ^(A p +Q - Ap _ Q J — U 



» 



++ 



(10) 



( v ++ 



+ V +_ 



-a!? } 



x(n) 
- A p , 



(A? 



(n) 



v +_ + v _+ - V _ j - o+ ~ V o J ~ u 



» 



» 



(n) 



++ 



4(A^ +0 + A^_ ) = 



,(«) 



(n) 



» 



4A B + (A p , ++ + A p , + _ + A p , _ + + Ap, __ + A p , 0Q - 



» 



(n) 



» 



P -0 



+ a!? } 



p o+ 



4- \ (n) ^ 



Here n ranges from 1 to oo, J\ = Jg, J n = J B for n > 2 and similarly for 
D n and K n , and A^ = 0. 

Letting n go to 00 in the equations above one obtains the bulk equations, 
which are only five, since in the bulk correlations like y 3 and y' 3 coincide. We 
have already solved the bulk equations in Ref. 13 and here we report the two 
coupled equations for the bulk order parameters 



Vi 



(B) 



[(v + -V-) + v (v*-v*)} /W 



V2 



{ 2 B) = [(V+ + V.)+ V (Vl + Vl) + 2 7 F + y_] jW 



(11) 



and the expressions for the correlation functions given as functions of the 
order parameters 



(B) 

y 3 

(B) 



[r}{Vl + V*) - 2 7 F + y_] /W 



vr = Wi - v>)} /w 

y { 5 B) = [r](Vl + Vl) + 2 7 F + y_] /W 
where 77 = exp(J B + K B ), 7 = exp(-J B + K B ), 



(12) 



V±=e 



(B) , (B) 



2(1 - y\ 



(B)s 



W = r]{Vl + Vl) + 2 7 V+V- + 2(V+ + V.) + 1 and y 



(B) 



(13) 



lim y 



(n) 



From equations (11) one can easily derive the equation for the second 
order critical temperature, which reads[13] 

1 5 



2C(7o - 1) 



C7o 

C(7o - 1) + 1 



(14) 



8 



where ( = e B cosh J B and 70 = 5 tanh J B . 

Turning to the whole set of equations (10) let us observe that, since is 
not possible to decouple different layers, some sort of approximation is in 
order; one can, for example, choose a maximum number n of layers and then 
solve numerically the equations for n = 1, 2, . . .n using the bulk solutions as 
boundary conditions, i.e. as if they were the solutions for layer n + 1. We 
choose the crudest approximation, n = 1, which allows us to obtain analytical 
results for the second order transition. The same assumption has been made 
in Refs. 5-8. 

The procedure we follow to determine the equations for the surface {n = 1) 
order parameters (which is analogue to that used in [13] for the bulk ones) 
can easily be shown to be equivalent to the natural iteration method (NIM) 
by Kikuchi[14], in the sense that the resulting equations, when solved by nu- 
merical iteration, give always thermodinamically stable solutions. Moreover, 
our procedure has the advantage of dealing with a considerably lower number 
of equations (only two, in the present case). 

From now on we will suppress the index which has become redundant; 
furthermore, instead of ^ we will write ( B \ because layer 2 now plays the 
role of the bulk. 

In order to determine the surface equations we express p v > = p p V and 
Pp — pp ^ as functions of p s = pp and of the bulk quantities; to this aim we 
consider a set of equations formed by the four equations in the set (10) for 
n = 1 which contain only elements of \ p > , the first two equations of that set 
for n = 2, and the three equations given by the first of the constraints (5) for 
p p '. This way we obtain a system which (recalling that A p / = hi/cy) can be 
solved for the elements of p p i , obtaining 



p p > Se = exp [5eJ B + (5e) 2 K B ] -j-p. 



6,e = +,0, -, 



(15) 



where 




5/6 



= 1 



(16) 
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and 

d s = J2 ex P i 6eJ B + (Se) 2 K B ] c e 5 = +, 0, -. (17) 

e 

The remaining five equations for n = 1, together with the normalization 
condition Tr(p p ) = 1 can now be used to obtain the elements of p p ; one finds 

p PSe = exp [5eJ s + (SefKs] -fsleG' 1 , 6,e = +, 0, -, (18) 

where 
and 



^"'^W \t ] -» = 1 (19) 



G = J2 [5e J 5 + (5e) 2 K s ] 7-5%- (20) 
Since the surface order parameters are related to p p by the equations 



Hi — Pp++ Pp — + Pp+o Pp—o 

2/2 = P P++ + P P + P P+0 + P P _ + 2p P+ - 

they must satisfy the following equations 

Vl = [e Js+Ks h 2 + -7^) + (7+ -7-)] /G 

y 2 = [e Js+Ks (~f 2 + + 7 2 ) + 2e- Js+Ks 7 +7 _ + ( 7+ + 7 _)] jG 



(21) 



(22) 



Once one has solved (22), using (8), (15) and (18) it is immediate to 
obtain explicit expressions for the n.n. two-points correlation functions. 

It is easy to check that, when y\ = 0, then y\ ' = too, that is, if 
the surface is paramagnetic, the bulk must also be paramagnetic. On the 
contrary, if the bulk is paramagnetic (y\ =0), the equation for y\ has 
always the paramagnetic solution y\ = 0, but can also have a ferromagnetic 
solution yi 7^ 0. It means that in the paramagnetic region of the bulk phase 
diagram one can have a surface transition line. The first order part of this 
line must be evaluated numerically, by comparison of the free energies of the 
two phases, while the second order one can be determined analytically by 
expanding (22) in powers of y\ (the order parameter y 2 must be expanded 
as well, and one must remember that y\ =0, because we are in the bulk 
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paramagnetic region of the phase diagram) around the solution yi = 0, and 
then requiring that the resulting equations are satisfied to the 3rd order in yi 
(going to the 5th order would give a set of equations for the tricritical point). 
One finds, for the surface second order transition 



2~ 8 e Ks -x 
Co 



3x + Vx 2 + 8 



-Kg _|_ 



X 



= e 



(23) 



where x = e Js - 2e~ Js > 0, c = 1 + 2V and c = 1 + 2Ve Ks cosh J B , with 



V = e" AB / 6 



(B) 
V2 



.2(1 -y ( 2 B) ) 



(24) 



3. Phase diagram and reentrant phenomenon 

In the present section we turn to the Blume-Capel model: thus we set 
K B = K s = 0. 

Following [5] we define the ratios R = Jb/Js and D = As /As and 
classify the possible phase diagrams at fixed R, D in the plane (d, r), where 
d = Ab/QJb and r = 1/6 Jb- As in [5] we obtain four main types of phase 
diagrams, which we report in Figs. 1-4; solid and dashed lines represent 
respectively second and first order transitions, and the thinner lines, when 
present, refer to the surface. B and S stand for bulk and surface respectively, 
while P and F stand for para- and ferromagnetic. 

The type- A diagram, reported in Fig. 1, is characterized by the presence 
of only ordinary transitions; the type-B diagram (Fig. 2) shows an ordinary 
transition at high temperatures, while in the low temperature region one has 
extraordinary and surface transitions: the three transition lines meet at the 
special point X; in the type-C diagram (Fig. 3) the situation is reversed, 
and finally, in the type-D diagram (Fig. 4) one has no ordinary transition 
but two non-intersecting, extraordinary and surface, transition lines. Type-B 
and type-C diagrams could be classified with even more detail, because in 
both cases the surface tricritical point can be either present or absent, but 
we do not enter into these point. 
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It is possible to distinguish between these four types by means of analytic 
conditions. First of all one must determine the bulk and surface critical 
temperatures at d = — oo, respectively denoted by r^ 00 "* and Tg°°' ) (the latter 
is determined under the assumption of paramagnetic bulk) , and the bulk and 
surface critical values of d at r = 0, respectively denoted by dg and d^\ 
and then compare these values, recalling that, if r^°°' ) > Tg.°°' ) or d^ > dg, 
the corresponding transition must be necessarily an ordinary one, because 
in the bulk ferromagnetic region of the phase diagram the surface is always 
ferromagnetic and surface transitions are forbidden. 

By means of (14) and (23) one obtains r^°° ) = (Slnf)" 1 « 0.822 and 
T ^ = (3RIH2)- 1 « 0.4811T 1 respectively; the condition r^ oo) = 4° o) 
gives the critical value 

R c = - 1 « 0.585, (25) 
In 2 v ' 

above which the d = — oo transition is ordinary (the corresponding MFA value 
is R c = 2/3). Furthermore, simple analytic considerations on the ground 
state yield d$ = 1/2 and d^ = from which the critical value 

D c = D C (R) = ^R (26) 

(equal to the corresponding MFA value), below which the zero temperature 
transition is ordinary, can be derived. Thus, the four main phase diagrams 
are characterized by the following conditions: 

Type- A diagram: R > R c and D < D c , 

Type-B diagram: R > R c and D > D c , 

Type-C diagram: R < R c and D < D c , 

Type-D diagram: R < R c and D > D c . 

Notice that type-B and type-C have not been found in the RG-based 
classification of Ref. 6, although a different renormalization group scheme[7,8] 
gives a phase diagram which is similar to our type B. 

To this fundamental classification we must add two new cases, which 
appear in the CVM treatment when D is just below D c and are illustrated 
in Figs. 5-6 for R above and just below R c , respectively. The new feature 
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of these diagrams is the presence of two or more special (i.e. multicritical 
points where surface and bulk transition lines meet) points. 

In Figs. 2,4 and 5 it is possible to observe a reentrant phenomenon in the 
ferromagnetic-paramagnetic (order-disorder) surface transition. In fact, for 
fixed R, there is a range of values of D for which one has such a phenomenon. 
The lower and upper limits depend on R, and the former is, of course, just 
below D c . In Figs. 7-10 we report the behavior of the order parameters and of 
the correlation functions for a set of values of the parameters corresponding 
to the dotted line in Fig. 4. As can be expected, the surface-bulk n.n. 
correlation functions are very small (even if non-zero) in the temperature 
range where the bulk is paramagnetic and the surface ferromagnetic. 

Finally, it is noteworthy that the surface reentrant phenomenon is not 
revealed by a MFA analysis: in Fig. 11 we compare the MFA phase diagram 
corresponding to Fig. 4 to that obtained in our approximation. 

4. The surface reentrance in the BEG model 

Let us now briefly consider the BEG model in the case £ > — 1, being 
£ = Kb/QJb (the case £ < — 1, where a staggered quadrupolar phase can 
occur, would require the introduction of two sublattices, with different order 
parameters to allow for the symmetry breaking between them) . The results 
of the previous section can be easily extended to this case; the only change 
that is needed in the classification is that one finds, for the critical value of 
the ratio of crystal fields, 

D < = l R TT^' < 27 > 

where V = Kb/Ks- Furthermore, since for the BEG model the bulk phase 
diagram can exhibit reentrance (see [13] and references therein) both for £ < 
and for £ > 3, it is interesting to know what happens in these cases for the 
semi-infinite model. 

For £ < one can easily find cases of reentrant ordinary transitions, 
where the bulk and the surface are simultaneously interested by the reentrant 
phenomenon, as in Fig. 12. 
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The most interesting situations occur for £ > 3. When D is slightly 
less than D c the surface reentrant transition already seen for the BC model, 
together with the ordinary reentrant transition, give rise to a double reentrant 
phenomenon, as shown in Fig. 13. Another double reentrant phenomenon can 
be obtained for D slightly greater than D c : in this case (Fig. 14) one has two 
successive reentrances on the surface transition line, the highest temperature 
transition being now second order. 

5. Final remarks 

We have investigated the semi-infinite spin-1 Ising models, devoting par- 
ticular attention to the BC model, in the framework of the pair approxima- 
tion of the CVM. We have classified the possible phase diagrams at fixed 
R, D, finding two new types of diagrams and showing that, at least in the 
present approximation, the first order surface transition do occur and ex- 
hibits reentrant phenomenon; finally, we have given qualitative conditions 
for this phenomenon and have discussed the extension of our results to the 
BEG model. 

The surface reentrant phenomenon we have found in the BC model is not 
revealed by a MFA analysis. It would be interesting to check our results by 
other high-precision methods. 

Our analysis has led to a classification scheme with four fundamental types 
of phase diagram, depending on the values of R and D. Two of them, namely 
type-B and type-C, have not been found in the approximations of Refs. 6 
and 10, while Refs. 7 and 8 present a phase diagram like our type-B. This 
point is worth of a further investigation, in order to establish whether the first 
order surface transition exists or it is only an effect due to the approximation 
method adopted. We conjecture that this transition is a real feature of the 
model, and give the following argument in support. 

Let us consider the BC model again: at T = 0, and for d > 1/2, the bulk 
ground state is paramagnetic, i.e. Sk = for all bulk sites k. Then all the 
interaction terms between bulk and surface vanish and the surface behaves 
like a BC model on an infinite square lattice, thus exhibiting a first order 
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transition (if D > D c ) at d = d s (indeed, this zero temperature transition 
is found also in [5]). Similarly, at a very low but finite temperature, only a 
small fraction of the bulk spins will be different from zero and the surface will 
be only weakly coupled to the bulk, this resulting in a slight change of the 
critical value of d, giving a smooth surface first order transition line in the 
(d, r) plane starting from the point (dg , 0). Work is in progress to confirm 
the conjecture and the results we have obtained. 
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Figure captions 



Fig. 1 : Phase diagram of the semi-infinite BC model for R > R c and 
D < D c . Solid and dashed lines denote second and first order 
transitions, respectively; X is the special point; B = bulk, S = 
surface, P = paramagnetic, F = ferromagnetic. 

Fig. 2 : Phase diagram for R = 1 and D = 1.8. Symbols as in Fig. 1. 

Fig. 3 : Phase diagram for R — 0.4 and D = 0.5. Symbols as in Fig. 1. 

Fig. 4 : Phase diagram for R — 0.4 and D = 0.8. Symbols as in Fig. 1. 

Fig. 5 : Phase diagram for R = 1 and D = 1.49. Symbols as in Fig. 1. 

Fig. 6 : Phase diagram for R = 0.56 and D = 0.8. Symbols as in Fig. 1. 

Fig. 7 : Behavior of the order parameters y\ (solid line) and 1/2 (dashed 
line) for the case corresponding to the dotted line in Fig. 4. 

Fig. 8 : Behavior of y 3 (solid line) and y' 3 (dashed line) in the case of Fig. 
7. 

Fig. 9 : Behavior of y 4 (solid line), y' 4 (dashed line) and y'l (dot-dashed 

line) in the case of Fig. 7. 
Fig. 10 : Behavior of y 5 (solid line) and y' 5 (dashed line) in the case of Fig. 

7. 

Fig. 11 : Comparison between the phase diagram of Fig. 4 (denoted by 
CVM) and the corresponding one obtained in MFA (denoted by 
MFA). 

Fig. 12 : Phase diagram of the BEG model for f = -0.5, R = 0.5, D = 0.4 
and r = 1. 

Fig. 13 : Phase diagram of the BEG model for f = 3.6, R = 0.5, D = 1.22 
and r = 1. 

Fig. 14 : Phase diagram of the BEG model for f = 3.6, R = 0.5, D = 1.26 
and T = 1. 
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